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ABSTRACT 

We present wide field photometry of the Galactic Globular Cluster M 22 in the B, 
V and I passbands for more than 186,000 stars. The study is complemented by the 
photometry in two narrowband filters centered on H a and the adjacent continuum, and 
by infrared J, H and K magnitudes derived from the 2 MASS survey for ~2000 stars. 
Profiting from this huge database, we completely characterized the evolved stellar 
sequences of the cluster by determining a variety of photometric parameters, including 
new photometric estimates of the mean metallicity, reddening and distance to the 
cluster. In particular, from our multi-wavelength analysis, we re-examined the long- 
standing metallicity spread problem in M 22. According to our dataset, we conclude 
that most of the observed width of the red giant branch must be due to differential 
reddening, which amounts to a maximum of AE(B — V) ~ 0.06, although the presence 
of a small metallicity spread cannot be completely ruled out. More specifically, the 
maximum metallicity spread allowed by our data is of the order of A[Fe/H]~ 0.1 -r 0.2 
dex, i.e., not much more than what allowed by the photometric errors. Finally, we 
identified most of the known variable stars and peculiar objects in our field of view. 
In particular, we find additional evidence supporting previous optical identifications 
of the central star of the Planetary Nebula IRAS 18333-2357, which is associated with 
M 22. 

Key words: globular clusters: individual: M 22 - planetary nebulae: individual: 
IRAS 18333-2357 



1 INTRODUCTION 

M 22 was one of the first Galactic Globular Clusters (GGC) 
to be discovered, in 1665, by Abraham Ihl e and also one 
of the first ones to be studied in detail JShaplevI fl93ot 
lArp fe Melbourndll95Sft . It s oon became the ta rget of a se- 
ries of studies (starting with lHesser et aI1ll977l) because of 
the large color spread of its red giant branch ( RGB) se- 
quence, similar to that observed in u Centauri JWoollevl 



* Based on observations made with the European Southern Ob- 
servatory telescopes, using the Wide Field Imager, as part of the 
observing program 65.L-0463. 

f This publication makes use of data products from the Two Mi- 
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of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation 
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1966). This suggested the possibility of a metallicity spread 
in M 22, as was demonstrated in the case of ui Cen a few 
years before dPickens fc Woollevlll967i) . 

However, while the presence of significant r eddening 
varia tions was excluded in the case of ui Cen llCannonl 
1980), some differe ntial reddening was found in the direc- 
tion of M 22 (see iRichter. Hilker. fc Richtlerl ll99St and 
references therein). Of course, the presence of differential 
reddening does not exclude the presence of some metal- 
licity spread, since the two effects could be both present 
and responsible for the observed width of the RGB. The 
phot ometric studies of M 22 by iPiotto fc Zoccalil (Il999l) 
and IRichter. Hilker. fc Richtlerl jl999T) were able to put 
upper limits to the amount of differential reddening of 
<7AKf y-n=0-05 and AE(B— V)=0.07, re spectively. Interest- 
ingly, IRichter. Hilker. fc Richtlerl lll999F) demonstrated that 
part of the spread observed in the Stromgren colours must 
be due to CH and CN variations, so that the eventual spread 
in heavy elements should be negligible. 
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Spectroscopic studies, on the other hand, gave con- 
troversial results as far as the heavy elements abun- 
dances are concerned, while a spread in the CH and 
CN abundances of RG B stars appears unquestionable 
jNorris fc Freeman! IT9831) 1 . For example, some studies re- 
ported on metallicity variations of 0.3-^0.5 dex in Ca 
and /or Fe, ofte n correlated with the CH and CN variations 
J^jetersoiJ^i^JPitachow^ 

Il222|; iLehnert. Bell, fc CoherJ Il99ll : iBrown fcW allcrstcin 
1992), while other studies found no significant varia- 
tion in the h e avy element con t ent jManduca fc Belj [l97S 
CoherJ Il98lt lOrattonl Il982l: lOratton fc Ortolanil Il989l 
Laird. Wilhelm. fc PetersorJll99ll) 

Concerning this a pparent contradiction, 

ILehnert. Bell, fc CoherJ l)l99lf) noted that, since the 
the estimated standard deviation of the abundance varia- 
tions (^0.2 dex) is close to the typical uncertainty of most 
spectroscopic analyses, it is very difficult to unequivocally 
demonstrate the presence of a metallicity spread. This is 
particularly true if one considers that most of the above 
studies are only based on a handful of stars 10). In the 
next years, thanks to the new generation of multi-object 
spectrographs, it will be possible to analyze large samples 
of stars in a homogeneous way, thus shedding more light on 
this issue. 

M 22 is a metal poor f[Fe/Hl~-1.62. lHarrisl and 
very bright cluster. Considering also its position on the sky, 
(1; b) = (9.89; -7.55), and its proximity to us (only 3.2 Kpc 
from the Sun 2 ) M 22 is certainly the ideal target for various 
studies, ra nging from the dynamical mod elling of dense stel- 
lar system ( Albrow, Dc Marchi, & Sahu 2002) to microlens- 
ing studies lSahu^^lTl200lh . On the other hand the char- 
acterization of its stellar content is still quite uncertain since 
it suffers, once again, from the presence of differential red- 
dening along the line of sight. Here we provide a complete 
and homogeneous photometric characterization of the stellar 
content of M 22. 

The paper is organized as follows. In Section [5] wc 
present the observations and data reduction procedures, 
compare our results with previous literature and derive op- 
tical and infrared mean ridge lines (MRL). In Section [3] wc 
deal with the differential reddening and metallicity spread 
issues. In Section 2] we derive the mean metallicity, redden- 
ing and distance along with other photometric parameters 
of M 22. In Section |5] we identify the known variables and 
peculiar objects, including the central star of the Planetary 
Nebula IRAS 18333-2357. In Section |S| we summarize our 
main results. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were obtained at the 2.2m ESO/MPI telescope 
at la Silla, Chile, using the Wide Field Imager (WFI), a mo- 
saic of eight 2048 x 4096 pixels CCDs. The instrument scale 
is 0'.'238 pix~\ giving a total field of view of 34' x 33'. A 
set of B, V and I images were secured during a single ob- 
serving run on 6-7 July 2000, with exposure times ranging 

1 Some marginal evidence for an over- abundance of s-process el- 
ements has also been reported by|(Jratton { 19821). 

2 Note that the Galactic Bulge is in the background of M 22. 



from 5 to 400 sec. We also secured a set of exposures in two 
narrowband filters, H a and 5ft, centered around the H a line 
(Ac ~6580 A) and on the adjacent continuum (A c ~6650 A). 
The un-calibrated H a photometry is briefly discussed in Sec- 
tion l5.ll The average seeing during the observations was ~1" 
full width at half maximum (FWHM). 

The raw images were corrected for bias and flat-field 
using specific IRAF 3 procedures, within the noao.mscred 
package. The photometric reduction was ca rried out usin g 
the DAOPhot II and ALLSTAR packages <StetsorJll987lb 
Stars were searched independently on each CCD of the WFI 
mosaic with a 3<r threshold and fitted with a second order 
spatially variable point spread function (PSF). 

We used standard IRAF routines to obtain aperture 
photometry for a sample of isolated stars at various posi- 
tions along the CCD. We derived the optimal radius for the 
aperture photometry by constructing the curve of growth 
of each star. We compared the aperture photometry with 
the DAOPHOT PSF-fitting photometry and we obtained 
an aperture correction of 0.00 for the V and I filters (with a 
typical error of 0.02 and 0.01 respectively) and 0.03 for the 
B filter (with a typical error of 0.02). The aperture correc- 
tions do not correlate with the position of the star on the 
CCD. 

Then we corrected our instrumental magnitudes consid- 
ering the extinction coefficients available for each filter from 
the ESO web page 4 and the airmass at the beginning of the 
observations. 

All observations were carried out under photometric 
conditions. The calibration to the standard Johnson- Cousins 
photometric system was obt ained using tw o standard fields 
(namely TPhe and PC 1323, lLandoltll992T) observed at dif- 
ferent airmasses during the night. The adopted calibrating 
equations are: 

B = & + 0.45 (b-v) -0.48 
V = v - 0.09 (6 - v) - 1.03 
I = i + 0.12 (v-i)-1.88 

where b, v and i are the corrected instrumental magni- 
tudes and B, V and I the corresponding magnitudes in the 
Johnson- Cousins photometric system. 

The resulting, calibrated color magnitude diagrams 
(CMDs) are displayed on Figures EJand^in the V,(B-V) 
and V,(V— I) planes, respectively. As can be seen, the popu- 
lation of M 22 dominates the CMD resulting from CCD #2, 
since the cluster center has been placed on that chip. The 
bulge and disk populations dominate instead the CMDs of 
the outer CCDs, where the contribution by M 22 tends to 
disappear (see, e.g., chips #4 and #5). A few bright stars 
with V<11.3, (B-V) = 1.8 and (V-I)=1.9 are saturated. 

Therefore, in the following sections, we will restrict our 
analysis to stars measured in CCD #2 only, since most of the 
cluster population lies in that chip, where the contamination 
by the disk and the bulge is less important and the CMD 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 

4 http:/ /www. eso.org/gen-fac/pubs/astclim/lasilla/index. html, 

see also: http:/ /www. Is. eso.org/lasilla/sciops/2p2/E2p2M/WFI/zeropoints/ 
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Figure 1. The V vs (B-V) CMDs obtained for each of the eight 
chips of the ESO Wide Field Imager mosaic. 
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Figure 2. The V vs (V-I) CMDs obtained for each of the eight 
chips of the ESO Wide Field Imager mosaic. 
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Figure 3. Th e M 22 photometry bvlPiotto fc Zoccalil il999l) (left 
panel) and bv lRosenberg et alJ J200(T fright panel) is superposed 
to the CMD obtained in this paper (light grey). 

appears cleaner. This way, we also avoid the propagation 
of the photometric zero-point differences between differ- 
ent CCDs, which could degrade the overall quality of the 
photometric catalog. 

Finally, in order to complement the multi-wavelenght 
study of M 22, we identified approximately 1840 stars in our 
catalogue having J, H and K 3 magnitudes measured by the 
2 MASS survey 5 . An example of the resulting Ks,(V— Kg) 
and K S ,(J-K S ) CMDs can be seen in Fig. 



2.1 Online Catalogue 

Although the following analysis is entirely based on stars 
belonging to CCD #2, several M 22 stars are still present in 
CCDs #1, #3, #6, #7 and #8 (Figures0and|3. Therefore, 
the full catalog presented in Tableland published electron- 
ically contains all the stars detected in these CCDs. 

The catalog contains the B, V, I calibrated magnitudes 
on columns #2, #4, #6, each followed by the formal daophot 
errors (<5B, <5V , SI). Only stars measured both in the V and I 
filters are tabulated while the B magnitude is provided only 
when available. A flag "0.000" can be found on columns #2 
and #3 if a star is not measured in the B filter. 

Star positions are provided both in the pixel coordinate 
system of each CCD (columns #7, #8) and in the equatorial 



Figure 4. The photometries of NGC 5272 (left panel), NGC 5897 
(middle panel) and the mean RGB ridge line of M 2 (right panel) 
are superposed on our M 22 CMD (light grey). See the text for 
references. 



(RA; Dec) coordinate system (columns #9, #10). In the first 
column we provide a sequential identifier for each star. 

The coordinates in the J2000.0 absolute astrometric sys- 
tem have been obtained with a procedure alr eady described 
in other papers (see e.g. iFerraro et alj|200ll) . The new as- 
trometric Guide Star Catalog (GSC II) recently released and 
now available on the web 6 , was used as reference. More than 
two thousand GSC II astrometric reference stars have been 
found in the field of view of each chip, allowing for an accu- 
rate absolute positioning of the image. In order to derive an 
astrometric solution for each WFI CCD, we used a program 
specifically developed at the Bologna Observatory (P. Mon- 
tegriffo et al 2004, in preparation). As a result of the entire 
procedure, r.m.s. residuals of ~ tf.'2, both in RA and Dec, 
were obtained. This value can be considered as a represen- 
tative uncertainty of the astrometric calibration procedure. 

The photometry of stars belonging to CCDs #1, #3, 
#6, #7 and #8 has been shifted to match the photometry of 
CCD #2, taking into account the different CCD responses. 
The correction was calculated by fitting the mean ridge lines 
calculated in H2.3l to the CMDs. The correction applied to 
each of the B, V, I magnitudes of the external chips is al- 
ways ^0.1 mag. However, in the CMD of chips #6, #7 and 
#8, only a handful of cluster stars is present in the main 



5 http:/ /www. ipac.caltech.edu/2mass 



http://www-gsss.stsci.edu/gsc/gsc2/GSC2home.htm 
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Table 1. CCD B, V, I online photometric catalog of M 22. Only a few entries are displayed to illustrate the catalog format and contents. 



Chip#2 


Star Id 


B 


<5B 


V 


<5V 


I 


<5I 




Y 


RA 


Dec 


1 


12.071 


0.010 


10.351 


0.010 


8.544 


0.124 


806.273 


1876.142 


279.11637069 


-23.92032070 


2 


10.854 


0.010 


10.483 


0.011 


9.938 


0.014 


1808.769 


4086.802 


279.04357491 


-23.77422799 


3 


12.632 


0.010 


10.835 


0.010 


8.925 


0.015 


1221.985 


1819.724 


279.08629030 


-23.92412172 


4 


12.707 


0.010 


10.885 


0.010 


8.936 


0.018 


1806.331 


1933.347 


279.04400418 


-23.91672994 


5 


12.793 


0.010 


10.941 


0.010 


8.979 


0.015 


52.285 


1092.889 


279.17106212 


-23.97201436 


6 


12.888 


0.010 


11.043 


0.010 


9.123 


0.015 


1406.258 


2075.037 


279.07292590 


-23.90726310 


7 


12.958 


0.010 


11.095 


0.010 


9.096 


0.018 


333.186 


757.506 


279.15076070 


-23.99417931 
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Figure 5. The photometries of NGC 5272 (left panel), NGC 5897 
(middle panel) and M 2 (right panel) are superposed on our M 22 
CMD (light grey). See the text for references. 

sequence and turn off regions and the calculated corrections 
are correspondingly less certain. 

2.2 Literature Comparisons 

In order to check our calibration, we compared our V,(V— I) 
photometry with two catal ogues previously publis hed by 
iPiotto fc Zoccalil Jl99Sf) and iRosenberg et al.l (boOOl) . 

Fig. [3] shows the resul ts of such a comparison . In par- 
ticular, the photometry bv lPiotto fc Zoccalil Jl999l) appears 
in good agreement with ours as far as the upper RGB is 
concerned (left panel of Fig. However, their horizontal 
branch (HB) appears redder and fainter than ours, by ap- 
proximately 0.18 dex in magnitude and 0.04 dex in color. A 
similar discrepancy can be observed for the lower RGB and 



the turn-off (TO) regions, pointing towards a possible resid- 
ual (V— I) color-term between t he two calibrations. On the 
other hand, the comparison with lRosenberg et al.l l)2000F) re- 
sults again in a discrepancy, but in the opposite sense (right 
panel of Fig. [3J . Both the RGB and HB appear in fact sys- 
tematically bluer and brighter than ours, by roughly 0.1 dex 
in color and 0.05 dex in magnitude. In this case, there ap- 
pears to be a residual zero-point difference between the two 
calibrations. 

To further investigate on this issue, we compared 
(Fig. our photometry with that of three globular clusters 
having similar metallicities (and thus similar RGB shapes) 
as M 22. More specifically, the left panel of Fig.^Jshows the 
comparison of our M 22 photometry with that of NGC 5272 
jFerraro et al.|l 9921). whi c h has a metallicity of {FefH.]zw=- 
1.66 in the lZinn fc West! (Il984l) scale. NGC 5272 has been 
corrected with the reddening and distance tabulated by 
iFerraro et alJ ill 9991) and with the corresponding values for 
M 22 (see Section 0}. The m iddle panel of Fig. shows the 
comparison with NGC 5897 (IFerraro et alJll9 97). which has 
[Fc/H]zh'=-1-68 and has been corrected for reddening and 
distance as above. Finally, the right panel of Fig. 0] shows 
the comparison of our M 22 photometry with the RG B mean 
ridge line of M 2, which has [Fe/Hlzw = -1.62 faarrisl 
Il996|) and was published bv lDa Costa fc Armandrofll (119901) 
already corrected for reddening and distance. We thus ap- 
plied t h e dist ance modulus and reddening of M 22 following 
Harris] (|l99£|), to be consisten t with the distance scale of 
Da Costa fc Armandrofll Jl99Cj) . 

As can be seen, all the three panels of Fig. [I] show an 
excellent match with the present photometry of M 22, thus 
dispelling any remaining doubt on the adopted absolute V 
and I calibration. 

Unfortunately, among the recent CCD studies of 
M 22, the only available B -band dataset is the one by 
iKaluznv fc Thompson! (1200 ill and, as the authors explicitely 
state in their paper, their absolute photometric calibration 
is not reliable. In Fig. [5] we c ompared our photo metry to 
that of NGC 5272, N GC 5897 iFerraro et al.lll999l) and M 2 
<Lee fc Carnevlll999l) in the V vs B-V plane. The CMDs 
of NGC 5272 and NGC 5897 have been corrected exactly 
as in Fig. 0] The M 2 CMD have been corrected for the 
appropriate reddening an d distance and the corresponding 
M 22 values tabulated by iHarria (1996). We find a reason- 
able agreement with each of the reference clusters and we 
conclude that also our B-band calibration can be considered 
reliable. In particular, a good match is obtained in the case 
of NGC 5897 and of the RGB of NGC 5272. The NGC 5272 
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HB is redder than the M 22 one even if the mean level ap- 
pears similar. We also obtain a reasonable agreement in the 
case of M 2, even if its CMD is somewhat bluer than the 
M 22 one. 

2.3 Optical and IR Mean Ridge Lines 

Profiting from this large photometric database, we have con- 
structed the RGB mean ridge lines for M 22, both in the 
optical and infrared colors. Since the optical photometry is 
deeper than the infrared 2 MASS photometry, we were able 
to reach down to the main sequence in the V,(V— I) and 
V,(B— V) planes, while we reach the base of the RGB in the 
K,(V-K S ) and K,(J-K 5 ) planes (see Fig.0. 

To derive the lines, we sliced the sequences in appropri- 
ate magnitude and/or color bins of variable size (depending 
on the number of points and on the shape of the sequence). 
The representative point of each bin has been computed as 
the 2 a clipped average. The set of representative points has 
then been fitted by analytical polynomials of variable degree 
until the best-fitting polynomial was found. The resulting 
mean ridge lines are tabulated in Tables [5] (optical) and [3] 
(optical-infrared) and they are overplotted on the respective 
CMDs, in Fig.0 

The morphology of the M 22 RGB in the in frared pass- 
bands was studied bv iDavidee fc Harrisl <ll996fl . These au- 
thors noted a pronounced discrepancy in the lower RGB be- 
tween their M 22 (K, J- K) fiducial line and that of M 13 
llDavid gefcH arris! Il995l) . a cluster of similar metallicity 
iHarriall996l) . We compared the shape of our (Ks, J-Ks) 
fiducial lin e (converted to the standard system using the 
relation of ICarpenterl|200ll) to the one of M 13 published 
bv ID avid^e^^Ha^ial^95ri and we find a reasonable agree- 
ment among the two ridge lines (see figure [(J. 



3 DIFFERENTIAL REDDENING AND 
METALLICITY SPREAD 

The most striking characteristic of the CMD of M 22 is the 
large color spread of the RGB, incompatible with measure- 
ment errors. As summarized in Section Q there is a long- 
standing debate about a possible metallicity spread in M 22, 
but the well established presence of som e differential red- 
dening llRichter. Hilker. fc RichtleiHl999f) in the direction of 
M 22 further complicates the analysis, since both mecha- 
nisms can contribute to widen the RGB in color and they 
are difficult to disentangle. 

The presence of differential reddening is easily demon- 
strated in Fig. |H| As can be seen, stars on the red and blue 
sides of the mean ridge lines derived in Section 12.31 occupy 
different spatial positions in the cluster. In particular, stars 
redder than the mean ridge line tend to populate preferen- 
tially the uppermost part of CCD #2 (i.e., North of the 
cluster center). The opposite is true for stars bluer than the 
mean ridge line. Thus, the northern part of the cluster must 
be more reddened than the southern part. The same kind of 
behaviour can be observed using the HB sequence, which is 
less sensitive to metallicity than the RGB, confirming that 
the dominant contribution to the color spread of both se- 
quences must be due to reddening variations. However, we 
still cannot exclude the presence of some (small) degree of 
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Figure 6. Our M 22 infrared mean ridge line (continuous 
line) is compared with the M 22 and M 13 normal points by 
iDayidge fc Harrisl ll996t) (open squares) and iDavidge fc Harris! 
11995) (filled triangles), respectively. We applied horizontal and 
vertical shifts to align M 13 to the M 22 normal points at the 
turn off and the point 0.05 mag redward of the turn off, just as 
done by Davidge fc Harris! il99ft) (see lower panel in their figure 
7). 

metallicity spread. From a careful inspection of Figure |H1 it 
is also evident that differential reddening is present on differ- 
ent scales, ranging from ~100 pix (see the "hole" at (X, Y) 
= (600, 1100) in the lower-left panel) to ~1000 pix (see the 
lower half of the lower- left panel). Therefore, a wide field 
study is also required in order to obtain mean properties 
that are really representative of M 22. 

Here, we profited from our multi-band photometry to 
put quantitative constraints to the amount of metallicity 
spread allowed by the color spread of the RGB. In particu- 
lar, we studied the RGB color distribution of M 22 in (B— V), 
(V-I) and (V-K s ), selecting stars with V^15 for (B-V) 
and (V— I), and stars with Kg ^12 for (V— Ks). This is neces- 
sary to avoid any contamination from the galactic bulge and 
to use only measurements with the highest possible signal- 
to-noise ratio. We then computed the color difference, at 
fixed magnitude, between each star and the corresponding 
mean ridge line (Section 12. 3H . The histograms of the result- 
ing differences are plotted in Fig.UUfor each color. A gaussian 
curve representing the measurement errors is overplotted on 
each histogram (dot-dashed curve). A gaussian fit to the 
actual color dispersion is also overplotted (solid curve) on 
each histogram. The latter gaussian shows of course a larger 
dispersion than the one representing the errors only, since 
it contains also the contribution of the intrinsic spread, i.e. 
the differential reddening plus the eventual metallicity dis- 
persion 7 . 

7 We point out the presence of a clump of stars with redder colors 
in Fig. [5] especially visible in (B— V). If the colors of these stars 
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Table 2. Optical M 22 mean ridge lines 



V 


(B-V) 


(v-i) 


V 


(B-V) 


(v-i) 


V 


(B-V) 


(V-I) 


V 


(B-V) 


(v-i) 


V 


(B-V) 


(V-I) 


11.25 


1 


752 


1 


928 


13.49 


1 


203 


1 


452 


15.68 


1.040 


1 


298 


17.47 





812 


1 


060 


18.57 


0.845 


1.117 


11.38 


1 


703 


1 


877 


13.60 


1 


190 


1 


444 


15.78 


1.034 


1 


293 


17.52 





809 


1 


059 


18.62 


0.851 


1.122 


11.49 


1 


659 


1 


830 


13.69 


1 


181 


1 


436 


15.87 


1.029 


1 


287 


17.57 





805 


1 


058 


18.67 


0.857 


1.127 


11.60 


1 


621 


1 


796 


13.78 


1 


171 


1 


429 


15.97 


1.024 


1 


283 


17.62 





803 


1 


057 


18.72 


0.862 


1.133 


11.69 


1 


593 


1 


763 


13.89 


1 


159 


1 


422 


16.07 


1.019 


1 


278 


17.67 





801 


1 


057 


18.78 


0.870 


1.140 


11.79 


1 


561 


1 


733 


13.99 


1 


149 


1 


415 


16.18 


1.013 


1 


273 


17.73 





800 


1 


057 


18.82 


0.874 


1.144 


11.89 


1 


528 


1 


705 


14.08 


1 


141 


1 


409 


16.27 


1.008 


1 


269 


17.77 





799 


1 


058 


18.86 


0.880 


1.149 


11.99 


i 
i 


^yo 


i 
i 


d i y 


14.18 


i 
i 


1O0 


i 




16.37 


1.002 


i 




17.82 


n 
u 


son 

OUU 


i 

i 


uoy 
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Table 3. Infrared and Optical-Infrared M 22 mean ridge lines 
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(J-K S ) 


6.814 


4.269 


1.020 


8.234 


3.723 
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We can then represent the observed color dispersion, 
a (M 1 -M 2 ) obs i as the sum of three terms, one due to pho- 
tometric errors, ^(a^-a^) 8 , the second to the differential 



are due to differential reddening, they could be tracing a denser 
interstellar matter region, with an E(B— V) which is ~0.06 mag 
higher than the average cluster reddening. 

8 Photometric errors are computed as the standard deviation of 
repeated mcasurmcnts of a star magnitude, available for the B, V, 
and I filters. In the case of the Kg filter, we used errors provided 
by the 2 MASS extraction algorithm. We excluded stars belonging 
to the inner 1' around the cluster center, where crowding effects 
are most severe. 



reddening, cr&EiMx-i/^)' an d the third to the intrinsic metal- 
licity spread, ct A [m/h] 

2 2 .2 2 

<J (M 1 -M 2 ) obs ~ a S(M 1 -M 2 ) + {T AE(M 1 -M 2 ) + a A[M / H] 

2 2 
= <Td(M 1 -M 2 ) + a (B-V) int 

Since the observed color spreads and the photometric errors 
are 

a (B-v) obs = 0.026; o s{B _ v) = 0.015 

<?(v-i) ohs = 0.030; °s(v-i) = 0.016 

a (.v-K s ) ob , = 0.055; cr S (v-K s ) = 0.030 



Wide Field Photometry of M 22 9 



ti p r 




Figure 7. The mean ridge lines listed in Table 121 and 151 are su- 
perposed to the CMDs in the optical (upper panels), infrared 
(lower-right panel) and optical-infrared planes (lower-left panel). 



Figure 8. The stars in the M 22 RGB are divided into two sam- 
ples, according to their position with respect to the mean ridge 
line (upper panels). The stars in the two sample have different 
spatial distributions (lower panels). 



we can derive the intrinsic color spreads obtaining: 



'(B-V) int 



J (V-K s ) int 



AE(B-V) 



+ a 



A[M/H] 



0.02 



AE(V-I) 



+ o 



A[m/h; 



= 0.03 



AE(V-K S ) 



A[M/H] 



= 0.05 



These values have to be considered as upper limits to 
the amount of differential reddening needed to explain the 
RGB width, if we assume zero metallicity spread. The above 
values are in reasonable agreement with pr evious determina- 
tions. In fact, <Jae(v-i) = 0.05, derived bv lPiotto fc Zoccalil 
( 1999) using the main sequence, compares reasonably 
with our a Aruv _ n = Qm, while AE(B-V) = (0.07;0.08) , de- 
rived respectively bv | Ric hter 1 Hilker. fc Richtlerl lll999l) and 
lAnthonv-Twarog. Twarog^^Crai?1^95l) from Stromgren 
photometry, compares very well with our AE(B — V)~ 3 • 

O"A£(B-V)=0.06. 

We also note that the intrinsic spread, in the three 
different colors, changes following the reddening laws (in 
the following we will always us e the reddening laws 
bv iDean. Warren fc Cousins! lll978l) and ISavaee fc Mathisl 
(1979)). In fact, assuming (7ab(s-v)~ 0.02 as above, we ob- 
tain o&e(v-i) and <^ae(v-k s ) that are virtually identical 
to the values derived above 



^AE(V-I) 
VAE(V-K s ) 



1.34 ■ & AE(B-V) 
2.72 • <Ta.e(B-V) 



0.03 
0.06 



Therefore, since the intrinsic width of the RGB scales from 
the optical planes to the optical-infrared plane exactly as 
expected from the reddening laws, we have to conclude that 
the room left for an intrinsic metallicity spread must be very 
small. 

To confirm this conclusion, we made use of two redden- 
ing free color indices 

E(B — V) 



(B-V) 



(B-V) 



E(V-I) 
E(B - V) 



(V-I) 



(V-K s ) 



Qbvk ^ " E{V-KsY 

and we plotted them in the top panel of Fig. 1101 The posi- 
tion of stars in the (Qbvi, Qbvk) plane should thus depend 
only on the intrinsic properties of the stars (i.e., tempera- 
ture and chemical composition). In particular, the spread 
around the mean locus should depend only on the eventual 
metallicity spread and on the photometric errors, since the 
Q-color indices are independent on reddening by definition. 

Therefore, we derived the mean ridge line of the locus 
in Fig. 1101 (thick line) with the same method used in Sec- 
tion 12.31 We measured the distance of each star from the 
mean ridge line and we constructed the histogram shown in 
the bottom panel of Fig. 1101 Two gaussians have been over- 
plotted: the one obtained by propagating the photometric 
errors to the (Qbvi, Qbvk) plane (dotted-dashed line) and 
the one that best-fits the distribution (solid line). As can be 
seen, the two gaussians are virtually identical. 
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Figure 9. The observed color difference distribution between 
stars in the RGB and the corresponding mean ridge line is plot- 
ted in B-V (lower panel), VT (middle panel) and V-Kg (upper 
panel). Two gaussian curves representing the fit to the observed 
distribution (continuous curve) and the distribution expected by 
photometric errors (dot-dashed curve) are also plotted on each 
panel. 



Summarizing, the metallicity spread contribution to the 
intrinsic width of the RGB must be smaller than our typical 
measurement errors, i.e., of the order of ~ 0.01 -f- 0.02 mag. 
Using equation 8 in ICarretta fc Braeaglial lll998l) and con- 
sidering (V-I)o, 9 =0.92 (see next section), we derive that a 
color spread of ~ 0.01 -r- 0.02 corresponds to a metallic- 
ity spread of approximately A[Fe/H]~ 0.1 -r- 0.2 dex, which 
is the maximum metallicity spread allowed by the present 
photometry. This spread is also of the order of the uncer- 
tainty of high- resolution abundance determinations for a sin- 
gle star (~0.15 dex). Therefore such a low metallicity spread, 
if present, could remain hidden into the instrumental errors 
even studying a large sample of high resolution spectra of 
M 22 stars. 



4 METALLICITY, REDDENING AND 
DISTANCE 

In this Section we characterize the photometric properties 
of the stellar population in M 22 by measuring the complete 
set of observables for the evolved sequences, with the aim of 
obtaining modern, CCD based estimates of the cluster mean 
metallicity, mean reddening and distance modulus. 

The Vzahb level was obtained from the comparison 
of the observed star distribution along the HB with syn- 
thetic HB of appropriate metallicity ([M/H] ~ —1.5, see 
below), following th e semi-empirical approach described by 
iFerraro et alJ (JT999) . This method is extremely robust since 
it allows the direct determination of Vzahb, instead of re- 
lying on the mean HB level as an indicator of the ZAHB 



Figure 10. RGB stars in M22 are plotted in the plane of the two 
reddening free color indices Qbvi vs Qbvk (upper panel). The 
continuos curve is the mean locus of stars in the plane. In the lower 
panel the observed distribution of the stars distances from the 
mean locus is plotted (histogram). The two gaussian curves rep- 
resent a fit to the observed distribution (continuos curve) and the 
distribution expected considering the photometric errors (dotted- 
dashed curve). 



level. We obtained 

Vzahb = 14.33 ± 0.05 
about 0.15 mag fainter than what reported bv lHarri^ 



1996). 
Harris 



of 



This i s easily understood in view of the fact that the 
( 1996) value is referred to the mean HB level instead 
the Zero Age Hori zontal branch level (see equation #2 in 
IFerraro et atlll999l) 

^osenbere et all feOOOD derived Vzahb =14.25 for 
M 22 and, considering the zero point difference with our V 
magnitudes (0.05, see section l2~2t . our Vzahb level results 
to be about 0.03 magnitude fainter than their value. 

Next, we proceeded to the simultaneous determina- 
tion of the mean metallicity and re ddening, us i ng a pro- 
cedure that was first introduced bv ISaraie dini (1994J). In 
its original descr iption, the m e thod allowed the derivation 
of [Fe/H] in the IZinn fc West! Jl984h scale, using V,(V-I) 
photometries and mean ridge lines of the RGB. Later, 
the p rocedure was calibrated to the I Carretta fc Grattonl 
l|l997l) scale bv ICarretta fc B ragaglial jl99ct) and exte nded 
to the B,(B— V) pla ne bv ISaraiedini fc Lavdenl lll997l) and 
IFerraro et all Jl99Sl) . 

Using the calibration of ICarretta fc Bragaj djajjl99S|). 
we assumed a mean HB level < V R r >=14.15 jHarrislll996ll 
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Figure 12. Differential (lower panel) and cumulative (upper 
panel) luminosity functions for stars selected on the RGB of 
M22. The arrows indicate the location of the RGB-bump at 
V=13.90±0.05. 



Figure 11. Infrared M22 CMD. The dashed line indicate the 
putative position of the RGB-Tip. 



and used the V,(V— I) mean ridge line derived in Section !^. 31 
obtaining [Fe/H] CG = -1.63 and E(B - V) = 0.36 (we also 
derived (V-I)q, o =0.92 and A Vi. 2 =2.19) 9 . Using the calibra- 
tion bv lFerraro et al ] Jl999l) . we adopted the Vzahb just 
derived and used the V,(B— V) mean ridge line, obtaining 
[Fe/H] CG = -1.73 and E(B - V) = 0.39 (we also derived 
(B-V) 0jS =0.73 and AVi. 2 =2.68) 

Since the two above reddening determinations are in 
good agreement with each other, we will adopt their average 
for the rest of this article: 

E(B-V) =0.38 ±0.02 

O ur determination of the reddening is identical to the 
one bv lRichter. Hilker. fc Richtlerl il999T) which is also the 
most recent determination. Other estimates of the redden- 
ing range from E(B— V)=0.32 to 0.42 a nd can be fou nd in 
[ Hesse J j 19761): lHarris fc Racmel (Il979l) : lHarrisI (I1996D and 
ICrockerl ([1988). These values are in broad agreement with 
our determination within the errors, if we take into account 
the presence of a strong differential reddening in direction 
of M 22 (see previous section) . 

The average metallicity of M 22 can also be derived from 
the optical-infrared CMDs V,(V— Ks) and Ks,(J— Kg) using 
the RGB intrinsic colors (V— Ks)o and (J— Ks)o, measured 



9 (V-I)o, 9 and (B-V)o, 9 are the dereddened RGB colors at the 
level of the HB, while AV1.2 is the difference in V between the 
HB and the RG B at a dereddened color 1.2 (i.e., (V- I)p=1.2 and 
(B-V)o=1.2, see lSaraiedinilll994tlFerraro et aill999l) . 



at different magnitude levels (s ee the definition of these pa- 
rameters inlFerr aro et alJl2000l). We used the calibrations of 
IValenti et al.| l|2004 1 in preparat ion!) to yield metallicities in 
the lCarretta fc Grat ton (1997) scale from the infrared pho- 
tometry in the 2 MASS photometric system. We corrected 
our optical-infrared CMD using E(B - V) = 0.38 and the 
(m— M)o derived below. This way we obtained different pho- 
tometric metallicity estimates for M 22, which produce an 
average value of: [Fe/H] CG = —1.68. 

Averaging together all the optical and infrared metal- 
licity determinations we obtained: 

[Fe/H] CG = -1.68 ±0.15 

This value is in agreement, within the errors, with the most 
recent spectroscopic and photomet ric works which provide 
mean metallicities of -1.55 , -1.48 l|Lehnert. Bell, fc Cohenl 



Il99ll ICarretta fc Grattonl 



jRichter. Hilker. fc Richtlerl 



1999) 



1997, respectively) and -1.62 



An estimate of the global metallicity can be now derived 
accord ing to the prescription of lSalaris. Ghiefii. fc Stranierol 
(1993): 



[M/H] = [Fe/H] ± log(0.638 ■ 10 [q/Fc] ± 0.362) = -1.47 

by assuming [a/Fe]=±0.30 l|Salaris fc Gassisil Il996l) . The 
distance modulus is derived from the comparison of the 
observed value of the ZAHB obtained above (Vzahb = 
14.33 ± 0.05) and the absolute lev el computed from the 

[Straniero 1 Ch ieffi 1 fc Limongil il997l^ models (see Eq 4 by 
Ferr^o^^d'^^fll ) . From this relation and the global metal- 
licity computed above ([M/H] = -1.47) we get My AHB = 
0.59. From this figure we derived an apparent distance mod- 
ulus 

(m-M)v = 13.74 ±0.2 
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and finally (using E(B — V)=0.38), an intrinsic distance mod- 
ulus (m-M)o=12.56 ± 0.2 which corresponds to a distance 
of ~ 3.2 Kpc. A conservative uncertainty of 0.2 mag is as- 
sumed. Our (m — M)y is abou t 0.15 mag fainter than the 
value reported bv lrlarrisl <ll99rj) . Therefore the two distance 
moduli are in agreement within the uncertainties , however 
a systematic difference between the iHarris! Jl996T) distance 
scale and distances derived by our sem i-empirical procedur e 
still remains and is further discussed in lFerraro et al.Nl999Tl . 

As a consistency check for the parameters derived until 
now, we derived the putative position of th e RGB-tip, in 
the K band. From the relations provided bv lFerraro et alJ 
(2000) and assuming the [Fe/H] cg just derived, we obtained 
as absolute magnitude of the RGB-tip: 



-5.93 



Then, using the (m-M)o and E(B— V) we derived the pre- 
dicted position for the RGB-Tip in apparent magnitude: 
K ti P = 677 

In Fig. El the K Sj (J-K s ) CMD of M 22 is plotted. 
Even if the presented CMD does not possess enough stars 
in the upper RGB for a safe determination of the RGB-Tip, 
it is nevertheless clear that the predicted RGB-tip is very 
similar to the magnitude at which the star counts drop to 
zero, Ks ~6.74. 

The last important feature that we were able to measure 
on the RGB of M 22 is the so-called RGB-bump. This fea- 
ture was one of t he first successful predictio ns of the stellar 
evolution theory jThomaslll967l : [lbenlll96Sf) . identified ob- 
servationally in 47 Tuc and subsequ e ntly in all the properly 
observed clus ters |King et alJll985 t iFusi Pe cci et al 199C ; 

m 



Ferraro et aTlll999t IZoccali et alJll99St iFerraro et allboOC : 
Cho fc Leell2002ft . both in the optical and infrared filters. 
Recently, the RGB-bu mp was also observed in a few satel- 
lites of the Milky Way jMonaco et alJl2002l : iBellazzini eta!] 
I2002L and references th erein) . 

As demonstrated bv lFusi Pecci et al] jl990F) . the change 
in slope of the integrated luminosity function is the safest 
way to identify the RGB-bump location, since it makes use of 
stars contained in several magnitude bins. We thus identified 
the RGB-bump of M 22 in the V band, as shown in Fig. E21 



V 



bump 



= 13.90 ± 0.05 



We also identified the RGB-bu mp in the Kg ban d, finding 
a value in good agreement with lCho fc Led <|2002T) . 

In Figure EH the differen ce of V bump and Vzahb , 
p versus [Fe/H] in the ICarretta fc Grattonl < llQ97f> 
scale, is plotted for the 42 clusters of iFerraro^t^ll jm)99l 
small empty squares) and for M 22 (big filled circle). As can 
b e seen, M 22 matche s quite well the empirical calibration 
of IFerraro etafl lll999D . 



• M22 (This Work) 

□ Ferraro et al. (1999) 




Figure 13. The magnitude of the M 22 RGB-bump (large filled 
circle) is compared with that of other globular clusters and with 
the empirical calibration of lFerraro et alJ ll999fl in the AVj^™ p 
vs [Fe/H] CG97 plane. 
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5 VARIABLE STARS AND PECULIAR 
OBJECTS 

M 22 is known to host various type of variab le stars, from the 
RRLyrae stars typical of globular clusters dWehlau fc Hoed 
Il97fet) to rather peculiar or rar e objects such as a 
type II cepheid, one sdB star dWehlau fc Hoed Il978t 
iKaluznv fc Thompson! 1200 J) and a Planetary Nebula (see 



Figure 14. The position of the variable stars recovered in our 
catalogue is plotted in the optical CMDs. Different classes of vari- 
ables have been plotted using different symbols. 
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Figure 15. (H a ; H a - 3?) CMD for M 22 stars lying within 150 
pixels from S* (small empty circles). S*, the optical counterpart 
of the PN in M 22, is represented as a large filled square, while 
stars in the HB phase are plotted as filled triangles. 



next section). A probable dw arf nova in outburst p hase has 
also been recently identified dAnderso n et al.ll200ot) as well 
as a population of X-ray sourceT dWebb^Gendn^fc Barretl 
2002). In spite of the modest number of Blue Straggler stars, 
M 22 con tains also a significant num ber of SX Phoenicis 
variables dKaluznv fc Thompsor]|200ll) . 

We cross-correlated our c atalogue with that of variable 
stars bv lClement et aP (1200111 which provides coordinates, 
classification and other useful information for all the known 
variable stars in M 22. In Fig. 1141 we used various symbols 
to show the position in the CMD of the various type of vari- 
ables successfully identified in our photometry. 51 variables 
over the 79 known in M 22 have been identified. Among 
them, 16 are RR Lyrae stars. Averaging the V magnitudes 
of the RR Lyrae observed at random phase in our photome- 
try we obtain Virs=14.17±0.25. This value for the HB mag- 
nitude is about 0.15 magnitude brighter than our estimate 
of Vzah b in q ualitative agreement to what expected from 
equation #2 bv lFerraro" et al.1 dl999h . 



5.1 The Planetary Nebula in M22 

Toge t her with Pal 6, NG C 6441 and M 15 (iJacobv et alJ 
Il997l: I Adams et ail 1 19841) . M22 is one of the few galactic 
globular clusters which are known to host a planetary neb- 
ula (PN). It was fir st discovered with the IRAS satellite by 
iGillett et ail dl986ft as a pointlike source (IRA S 18333-2357) 
and then identified as a Planetary Nebula dGillett et alJ 
1989) . The central star of the PN has also been identified 
(Gille tt et alJ |l989). as a blue star (the "southern compo- 
nent", hereafter S*) belonging to a pair which lies only ~ 2" 
away from the infrared source. 

Some of the physical properties of the PN an d its cen- 
tral star were derived bv ICohen fc Gilletti lll989l) . Consid- 



ering the expansion velocity (11 km/s) of the nebula and 
assuming that its size corresponds to the distance to S*, 
ICohen fc Gilletti dl989l) concluded that the age of the PN 
appears to be only ~ 6000 yr. This short time scale implies 
that the central star should still be quite bright, with a lu- 
minosity comparable to that of the RGB-tip. This poses an 
apparent problem, since S* has a magnitude (and even col- 
ors) comparable to that of a slightly evolved HB star, as can 
be seen in Fig. 1141 

To further investigate the photometric properties of S*, 
we constructed the instrumental (non-calibrated) H a ,(H a — 
K) CMD for stars contained in a circle with a radius of 
150 pixels, centered on S* (Fig. 1151 . In this diagram, RGB 
stars occupy a vertical sequence while genuine HB stars 
(filled triangles) , which have a more pronounced H a absorp- 
tion, tend to have redder (H a — Jft) colors. In the same dia- 
gram, S* (large filled square) shows a moderate H a excess, 
with (H a — 5ft) slightly bluer than the normal RGB stars and 
very different from the HB stars. This is in agreement with 
the substantial amo unt of hydrogen emission disp layed by 
the spectrum of S* dHarrington fc PaltoglovJll993Tl . consis- 
tent with what expected if S* is the central star of a new- 
born PN. 

We thus need to explain the unusual position of S* in 
the V,(B— V) diagram. An enlargement of the CMD around 
S* is shown in Fig. 1161 (upper panel), where the theoretical 
isochrone of a post-AGB st ar with t~12.5 Gyr and Z=0.0004 
has been also overplotted dBertelli et aljlT994l) . Clearly, the 
position of S* (large filled square) is not compatible with 
the isochrone. However, S* should be surrounded by nebular 
dust, whi ch causes an excess of reddening of about AE(B- 
V)~0.14 dCudworthl Il990l: lHarrington fc Paltogloul Il993ft . 
After applying the corresponding correction, the position 
of S* is reconciled with the isochrone (large filled circle in 
Fig. 1161 upper panel). 

The theoretical (log T e ff, log L/Lq) plane can also 
provide useful insight into the evolutionary state of S* 
and its compatibility with a post AGB-star of about 
6000yr. Spectroscopic studies deri ve for S* a tempera- 
ture of 50,000i;T ef f ^75,00 K dCohen fc Gilletd Il989l : 
lHarrington fc Paltoglovj Il993l) . In this range of temper- 
atures, the bolometric correction in the V band, BCv, 
varies between -4.08 and -5.30, using a blackbody model 
dOriglia fc Leithererl Eoool) . If we assume T e// =60,000 K 
(hence BCy=-4.62), AE(B-V)~0.14 and (m-M) v =13.75, 
then S* occupies the position marked by a large filled circle 
in the lower panel of Fig. 1161 We also overplotted the theo- 
retic al post-AGB evolutionary t rack for a 1 Mq and Z=0.001 
star dVassiliadis fc Woodlll994l) . Tiny empty dots mark the 
positions occupied by stars of 5500, 6900 and 8300 yr. As can 
be seen, S* closely matches the track at a luminosity com- 
parable to that of the RGB-tip and a temperature which is 
fully compatible with the central star of a PN of ~6000 yr. 

We thus have been able to add supporting evidence to 
the identification of S* as the central star of the planetary 
nebula IRAS 18333-2357 in M 22. 



6 SUMMARY AND CONCLUSIONS 

We presented a wide field (33' x 34' ) , multi-band photometry 
(B, V, I, H Q and the adjacent continuum, K) of the globular 
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our photometric measurements, most of the intrinsic width 
of the RGB must be due to differential reddening, while 
the maximum metallicity spread allowed by our data is 
A[Fc/H]~ 0.1 ~t~ 0.2 dex, i.e., compatible with the photo- 
metric errors. 

We finally identified most of the variable stars and pe- 
culiar objects known in the observed field of M 22. In partic- 
ular, we provided additional evidence supporting the opti- 
cal identification of the central star in the planetary nebula 
IRAS 18333-2357, one of the few identified in a GGC up to 
now. 
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Figure 16. Upper panel: zoomed CMD in the HB region. The 
position of S* is marked as a filled square. The filled circle repre- 
sents the putative position of S* if a reddening excess of AE(B- 
V)~0.14 is assumed. The continuous curve is a post-AGB theoret- 
ical isochrone with t~12.5 Gyr and Z=0.0004. Lower panel: The 
position of S* in the absolute plane, once applied the bolomet- 
ric correction, is plotted (filled circle). The continuous curve is a 
theoretical post-AGB evolutionary track for a 1 Mq and Z=0.001 
star. Empy circles mark evolutionary steps between 5000 and 8300 
years. 



cluster M 22. For the H Q and K filters we presented only 
instrumental magnitudes in a tiny area around the PN in 
M 22. We provided the astrometrically calibrated catalog 
containing B, V and I calibrated magnitudes for ~140,000 
stars covering an area of ~24'x33'. About 2000 stars in 
our catalog have been also measured in the near infrared 
by the 2 MASS survey, providing calibrated J, H and Ks 
magnitudes in addition to our five optical bands. 

We use this catalog to characterize the evolved stellar 
sequences in M 22, especially the RGB, by deriving mean 
ridge lines in the optical and infrared colors, and by mea- 
suring the V magnitude of the RGB-bump, which appears 
in good agreem ent with the most recent calibration ver- 
sus metallicity l)Ferraro et al.lfl999|l . We also derived the 
mean metallicity, reddening and distance moduli, in a self- 
consistent way. The derived values agree well with previous 
determinations. 

Profiting from our multi-band catalogue, we re- 
examined the problem of the metallicity spread in M 22. 
We demonstrated the pr esence of differential reddening, 
confirming earlier findings jRichter. Hilker. fe Richtlerll 999. 
and references therein). We also conclude that, according to 
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